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We discuss the sensitivity to anisotropies of stochastic gravitational- wave backgrounds (GWBs) 
observed via space-based interferometer. In addition to the unresolved galactic binaries as the most 
promising GWB source of the planned Laser Interferometer Space Antenna (LISA), the extragalactic 
sources for GWBs might be detected in the future space missions. The anisotropies of the GWBs 
thus play a crucial role to discriminate various components of the GWBs. We study general features 
... of antenna pattern sensitivity to the anisotropies of GWBs beyond the low-frequency approximation. 

\ We show that the sensitivity of space-based interferometer to GWBs is severely restricted by the data 

. combinations and the symmetries of the detector configuration. The spherical harmonic analysis 

' of the antenna pattern functions reveals that the angular power of the detector response increases 

with frequency and the detectable multipole moments with efi'ective sensitivity hss ~ 10 ^"Hz 
^ , may reach I ~ 8-10 at / ~ /, = 10 mHz in the case of the single LISA detector. However, the cross 

■ correlation of optimal interferometric variables is blind to the monopole {I — 0) intensity anisotropy, 

' and also to the dipole {£ = 1) in some case, irrespective of the frequency band. Besides, all the 

self-correlated signals are shown to be blind to the odd multipole moments {£ = odd), independently 
of the frequency band. 
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I. INTRODUCTION 

' Space-based gravitational wave detectors retain many possibilities of providing access to new gravitational-wave 
I sources that are not covered by ground-based gravitational-wave detectors. The Laser Interferometer Space Antenna 
(LISA) is such a planned gravitational-wave observatory aimed at detecting and studying low-frequency gravita- 
^ ' tional waves in the band 0.1 mHz ^ O.lHz. The constellation of the LISA and the next generation detectors, e.g. 
(DJO! DECIGO/BBO 0,11, will consist of three spacecrafts keeping a triangle configuration. 

' Compared to ground-based detectors, a space-based gravitational-wave detector is characterized by many different 
. ^ ! features. For instance, the LISA introduces complications unknown to ground-based detectors, such as the coinplex 
' signal and noise transfer functions. The complications block the analytical characterization of the detector j 
J-j [ and in particular, the response of an interferometer becomes complicated for gravitational waves shorter than the arm 
■ - - I length of the detector. Only in the low-frequency limit, the response of the LISA detector is very simplified. The 
three arms of the LISA function like a pair of two-arm detectors, and it is well known that the pair is equivalent in 
the low frequency limit to two 90° interferometers which are rotated by 7r/4 with respect to each other (e.g. |^). 

An important ingredient of a space-based detector is analysis of time-delayed combinations of data streams, which 
provide laser-noise-free interferometric variable. The technique to synthesize data streams is known as time delay 
interferometry (TDI) [Hill [11 El. Several TDI signals, such as Michelson-like and Sagnac-like signals which are 
free from the laser frequency noise, will have different responses to secondary phase noise sources and to incoming 
gravitational waves. Starting with the original TDI observables for stationary-array combinations, the TDI observables 
have been developed until recently (see [Tj and references therein). 

Space-based gravitational-wave detectors could be the most suitable devices to study and search for stochastic 
gravitational-wave signals. Examples of stochast ic g ravitational waves are those produced by large populations of 
Galactic |0, [l6l ITtI I and extra-galactic binaries |l8l |20| and a primordial gravitational-wave background pro- 
duced by several cosmological mechanisms (see |21| for a review). Stochastic gravitational waves are expected to be 
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anisotropic, and an important issue is to identify unambiguously the anisotropy to get insights into the origin and 
underlying physics of them. 

A method to explore an anisotropy of gravi tational-wave background has been recently proposed based on the time 
modulation of the single data stream 22, "2^, "2?, "2^, and/or the two data streams 24, 27, 28], which allow us to 
extract the individual coefficients of multipolc moments related to a distribution of gravitational waves on the sky. 
Hence, provided all the coefficients of m ultip olc moments obscrvationally, one can, in principle, make the sky map of 
the gravitational-wave backgrounds p3 . l25l |. It was demonstrated in the low frequency limit that the LISA is blind 
to the whole set of odd multipolc moments and sensitive only to monopole = 0), quadrupole {£ = 2) and octupole 
{£ = 4) anisotropy p^. |2^ . Actually, the multipolc moments £ = 2,4 and their m th harmonics of the galactic 
distribution of binaries would be observable with sufficiently high signal-to- noise ratios, except for some multipolc 
harmonics |27j . The restricted sensitivity to the multipolc moments is an immediate outcome of the low-frequency 
approximation, but what is the underlying physics that determines the limitation? As discussed in this paper, it is 
intimately associated with the geometric properties of the spacecraft configuration. 

So far most of the works aimed at probing the anisotropy of gravitational-wave background by means of space-based 
interferometers have been restricted to the low-frequency approximation. One reason is that a confusion gravitational- 
wave background formed by the superposition of many Galactic binaries comes in the low-frequency band of the LISA 
as the dominant source [l5|. However there are many expected sources of gravitational- wave backgrounds that 
spread outside the low- frequency region [isL Il9l l20l l29l l30| . Thus it is an interesting problem to create maps of the 
gravitational-wave background in a very wide range of frequency. For that purpose we need to know the general 
response and properties of space-based interferometer over a wide range of frequency. 

In this paper we are interesting in general features of response function for space-based detectors. The sensitivity of 
space-based detectors to multipole moments of a gravitational-wave distribution is in general restricted by symmetries 
of a response function independently of frequency. For example, symmetries of a response function tell us that a self- 
correlated data is blind to the odd multipole moments of anisotropy irrespective of frequency band (Sec. IIV Cp . Other 
interesting features independent of frequency can be also derived based on symmetries of a detector's response and 
geometric configuration of the spacecrafts. 

The paper is organized as follows. After briefly reviewing the detection method of an anisotropy by the correlation 
analysis in the next section, detector response functions for space-based interferometers are given in Sec. lIIII In Sec. lIVI 
we develop spherical harmonic analysis of antenna pattern functions and derive various fundamental properties of 
multipole moments. Based on those fundamental properties, in Sec.0 we examine the directional sensitivity of space 
interferometer. Angular power and effective sensitivity curves are discussed there, specifically focusing on the LISA 
detector. Section IVTl concludes the paper with a brief summary. Below the speed of light is set equal to unity (c = 1). 



II. DETECTION OF ANISOTROPY THROUGH THE CORRELATION ANALYSIS 

We begin by discussing how one can prove the anisotropy of gravitational- wave background based on the correlation 
analysis. A stochastic background of gravitational waves can be expressed as a random superposition of plane waves 
propagating along SI = — n direction with surfaces of constant phase ^(x) = t -I- n • x. Then the metric perturbation 
in transverse-traceless gauge h is expressed as: 

h(t,x) ^ J2 r df [dn e'^^fihA{f,n) e^in), (1) 

where J d^t denotes an integral over the sphere and hA{—f) = h\{f) are the Fourier amplitudes of the gravitational 
waves for each polarization mode. The Fourier amplitudes Ha is assumed to be characterized by the Gaussian random 
process: 

l^hAif,^)) = 0, 

(h*A{f,n)hA'{f,n')) = U{f ~ r f^^^^'h AA'SniUl^). (2) 

where 5*^(1/1, ^2) is the power spectral density of gravitational waves. The polarization tensors e"^(ll) appearing in 
equation ^ may be explicitly given as follows: 

e+(r2) = u(g) u - V (g) V, e^(r2) = u (g) V + V (g) u, (3) 

where the unit vectors u, v are expressed in an ecliptic coordinate as: 

u = cos cos -I- cos sin — sin 6'£;z, 
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V = sin(/)£;x — cos^ey, 

n = siuOe cos (j)EX + sinOEsincpEy + cos6e2 ^ —ft. (4) 

The detection of a gravitational- wave background is achieved through the correlation analysis of two data streams. 
The output signal for the detector / denoted by si(t) is described by a sum of the gravitational- wave signal hi{t) and 
the detector noise nj{t): 

si{t) = hi{t)+ni{t). 

We assume that the noise ni{t) is treated as a Gaussian random process with zero mean and spectral density Sn{f)- 

(niif)) = 0, 
{n}{f)njin) = \s{f-nSjjS^i\f\). 

Here, we further assume that the noise correlation between the two independent detectors is neglected ^. 

On the other hand, in addition to the information of a gravitational- wave background, the output signal hj(t) 
contains the time variation of the detector response caused by the detector motion. For example, the rotation of the 
Earth sweeps the ground-based interferometer across the sky. As for the space interferometer, LISA, the antenna 
pattern sweeps over the sky as the LISA constellation orbits around the sun with period of one sidereal year. These 
effects induce the signal modulation, which can be used to extract the information of anisotropy of gravitational-wave 
backgrounds. 

According to Ref.j25, we introduce two time-scales, AT and Tq: the light travel time AT between the two detectors 
(space crafts) and the period of the detector motion To. Since AT ^ To, it is possible to choose the averaging time 
scale r as AT <C t ^ To appropriately. Then, one can safely employ the correlation analysis between two detectors 
as a function of time averaged over the period r. Keeping this situation in mind, the output signal hi{t) may be 
written as 

/ + 00 I- ^ 

where the colon denotes the double contraction, i.e., D: e Dije^^ 31]. The quantity D is detector's response 
function, whose explicit expression will be presented in next section. Note that the response function depends on time 
due to the detector motion. 

Provided the two output data sets, the correlation analysis is examined depending on the strategy of data analysis, 
i.e., self-correlation analysis only using the single data stream or cross-correlation analysis using the two independent 
data stream: 

C{t) = {sij^t)sj{t)) = {hj{t)hj{t)) + {ni(t)nj{t)) 

= IZ^I ^S,{\f\,n)Tfj{f,n;t)+6jj |^|5„(|/|), 

where Tfj is the antenna pattern function defined in an ecliptic coordinate, which is expressed in terms of detector's 
response function and an optimal filter: 

Ff{nj-t)=T>i{n,f-t):e^{n) (6) 

with Q{f) being the Fourier transform of the optimal filter. Note that the phase factor g^'^'^ f ^-i^i-^j)/ ^ arises due 
to the differences between the arrival time of the signal at each detector. The above expression implies that the time 
series data C{t) as observable is given by the all-sky integral of the spectral density or luminosity distribution 
of gravitational waves convolving with the antenna pattern function. To see this more clearly, for the moment, we 



In the case of the space interferometer, while the various data streams can be constructed combining the signals extracted from respective 
space crafts, most of them are dominated by a correlated noise. Thus, the optimal data combinations which cancel the correlated noise 
are required to work with the correlation analysis. 
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neglect the noise contribution and set the optimal filter as Q{f) = 1. Keeping the assumption AT <^ t <^ Tq, the 
detector output C{t) is written as 

= f^^J§s,{\f\M):Ffj{f,t,n). (7) 

We then decompose the antenna pattern function J-fj and the luminosity distribution into spherical harmonics in an 
ecliptic coordinate, i.e., sky- fixed frame. We have 

SM n) = Y^ipLifW YLm, ^FAf, t) = Y.af^{f,t) yu^). (8) 

£.rn Irn 

Substituting ((HJ into becomes 

^(^'/) = i E [pfm{f)Y^fra{f, t) ■ (9) 
hn 

Note that the time dependent multipole coefficient af^ appears due to the detector motion, which can be eliminated 
by further employing the harmonic expansion in detector's rest frame {0,(j)). We denote the multipole coefficients 
of the antenna pattern in detector's rest frame by agm (see Eq.|5B]). The transformation between the detector rest 
frame and the sky-fixed frame is described by a rotation matrix by the Eulcr angles {ip, ip), whose explicit relation 
is expressed in terms of the Wigner D matrices ■ 

e 

aLif,t) = e-'"^dLWe-""^ (10) 

n=-l 

Here the Euler rotation is defined to perform a sequence of rotation, starting with a rotation by tjj about the original 
z axis, followed by rotation by about the original y axis, and ending with a rotation by if about the original z axis. 
Note that the Euler rotation conserves the multipole moment £, but mixes m-th harmonics. 

For illustration, let us envisage the orbital motion of the LISA constellation. The LISA orbital motion can be 
expressed by "0 = —Lut, "i? = — 7r/3, ip = wt, where lu = 2tt/To is EISA's orbital frequency (Tq = 1 sidereal year). Since 
the antenna pattern function is periodic in time due to the orbital motion, one can naturally perform the Fourier 
transformation of the detector output by [2^ |23 | : 

J^o Jo 

Using the relation IjlOl) . we finally obtain: 

-. oo £—k 

= ^T.ll ipfmifF W ae,im+k){f), (11) 

i=0 m=-l 

for > 0. The above equation shows how the detector output depends on the multipole coefficients aim in detector's 
rest frame for a given luminosity distribution of a gravitational- wave background, vFrXf)- Given the output data 
Ck{f) experimentally, the task is to solve the linear system Hll|l with respect to pf^if) if we know the antenna 
pattern function. As discussed by Cornish 24], this deconvolution problem is typically either over-constrained or 
under-constrained depending on the antenna pattern. In this sense, the understanding of the general properties of 
antenna pattern functions is primarily important and would shed light on the deconvolution problem. It might be 
further helpful to characterize the directional sensitivity to the sky map of the gravitational-wave background. The 
detailed investigation for the spherical harmonic analysis of antenna pattern will be presented in Sec. IIVI Before 
developing the analysis, we briefly review the detector response functions for space interferometers. 



III. DETECTOR RESPONSE FUNCTION FOR SPACE INTERFEROMETER 

In this section, according to the treatment based on the coordinate- free approach in Ref. [^l33j. we derive various 
types of detector response functions for space interferometer, which will be used to analyze the sensitivity to an 
anisotropy of gravitational-wave background. 
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A. One-arm detector tensor 

Following the Doppler tracking calculations described in Ref . ,31j , the optical-path length between spacecraft i and 
spacecraft j is formally written as 



(12) 



where g^j^ is spacetime metric. According to Ref. 6J, the optical-path variation in presence of the gravitational waves 
is given by 



6e^J{U) = £,j{U) [ df f dn D(/,t„n):h(/, f2)e* 2. 

J —00 J 



(13) 



where the one-arm detector tensor D and the transfer function T are 

^1 



T(/,tj,n) = sine 



2^ 



[l + n-r,,(i,)] 



exp 



*;r;^[l + n-ry(ii)] 



(14) 
(15) 



where rij{ti) is the unit vector pointing from the space craft i at the time of emission ti to the space craft j at the 
time of reception tj, i.e., rij{ti) = {x.j{tj) — Xi{ti)}/lij{ti). The function sinc(a::) is defined by sinc(a;) = s'mx/x and 
the variable /y = [2Tr£ij{tj)]~^ means the characteristic transfer frequency. 

B. Detector response function 

The calculation of the one-arm detector tensor can be applied to derive the response function for a space interfer- 
ometer via Doppler tracking method. The constellation of the planned space interferometer, LISA and also the next 
generation detectors DECIGO/BBO constitutes three space crafts and each of them is separated in an equal arm 
length (Fig. Note cautiously that the detector arm length varies in time, mainly due to the intrinsic variation by 
the Keplerian motion of three spacecrafts and the tidal variation caused by the gravitational force of Solar System 
planets [s^ . Issf . The caveats concerning these effects have been already mentioned "s^ and their influences were 
recently investigated. As long as the low-frequency gravitational waves with frequencies comparable or lower than 
the characteristic frequency fij are concerned @, l35l l3^ , the so-called rigid adiabatic approximation |35j| , in which the 
three space crafts rigidly orbit around the sun under keeping their configuration, really works. 

Keeping these remarks in mind, we adopt the rigid adiabatic approximation to give an analytic expression for 
response functions. For the sake of the brevity, we work with the static and the equal-arm limit of the detector 
response. In this case, one writes L = lij and /* = fij = l/(27rL). Thus, the interferometric signals combining with 
the six data streams can be generally expressed as function of 



/ 



(16) 



Specifically, for LISA detector, the arm length is L 5 x lO^km, yielding ~ lOmHz. 

Based on the configuration in Fig. ^ a signal of Michelson interferometers extracted from the space craft 1 is P.l33j: 



= —[Siuit - 2L) + S£2i{t - L) - Shsit ^ 2L) ~ Shiit - L) 



df / dfl DM,(l^,/):h(/,f2) e 



(17) 



where Dm^ (J^, /) is the detector tensor. The explicit form of the detector tensor is given by 



DM,(fl,/) = -{(a®a)rM(a.O,/)-(c®c)TM(-c.O,/)} 



e-^f \ sine ( ^^^""'"^ 1 e"3/(i+-n) ^ ^^^^^ 



/(I + U ■ ] g-|/(-l+u.f2) 



(18) 
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The directional unit vectors for the three space crafts are denoted by a, b,c (Fig. QJ. Note that the above expression 
possesses the symmetry, i.e., DM(f2,— /) = Dm(J1,/)*. 

Unfortunately, the simple Michelson interferometry with unequal armlengths does not cancel the laser frequency 
noise, which is thought to be one of the most dominant sources in the instrumental noises. Thus, the Michelson signal 
might not be a viable interferometric variable. Instead, a number of so-called TDI variables that cancel the laser 
frequency noise even when the armlengths arc unequal have been proposed . These signals are built by combining 
time-delayed Michelson signals so as to reduce the overall laser frequency noise down to a level of other noises. A 
particular example of a TDI variable is the X signal: 

^ Isiuit - 2L) - 5h2{t - 4L) + 5t2i{t ~ L) - «2i(i - 3L) 



AL 



-Sei3{t - 2L) + d£i3{t - AL) - Se^iit -L) + Sisiit - 3L) 



(19) 



This signal is expressed by a superposition of the Michelson signal, Sx-^{t) — ^[sf{t) — sf{t — 2L)]. Thus, the detector 

tensor for the interferometer variable is Dx = ^(1 — e~^*'^)DM. 

Other useful combination comes from comparing the phase of signals that are sent clockwise and counter-clockwise 
around the triangle. Such combination is named as Sagnac signal. The Sagnac signal extracted from the space craft 
1 is 

^ \sii3{t ~ 3L) + Si32{t - 2L) + 6£2i{t - L) - Si^it - 3L) ~ «23(i - 2L) - S^iit - L) 



3LI 

r df I dn Ds, {nj): h(/, e2-/«(-^) , 

J — oo J 



(20) 



where the detector tensor Dgj^ is expressed as 

1 



Ds,(f^,/) 

Uf) 
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{(a®a)ra(/) + (b®b)rb(/) + (c®c)rc(/)} ; 



e-^'/(-2+a-i^)sinc 



^(l + a-O) 



g-ii/(3+(a-c)-n) 



|(l + b.n) 



|(l-b.n) 



-(l-a-n) 

2 ^ ' 



(1 + c-O) 



/ 



(1 - c • n) 



(21) 



The three Sagnac signals extracted from the space crafts 1, 2 and 3 are often quoted as a, /?, and 7 in the literature (e. 
[ll|). Combining these variables, a set of optimal data combinations free from the noise correlations is constructed fl 
(see also Q): 



= "^('^si -2/1S2 +^83), 



+ /iso + hs„). 



(22) 



It is worthwhile to note that in the low- frequency limit / <C 1, the detector tensor for the Michelson, the X and 
the Sagnac signal can be simply expressed as 



3 i f 

ifB^,^ = Dx, = - Ds, = y (a®a-c®c) 



Oif) (23) 
Using the above expression, the detector tensors for optimal combinations A, E and T respectively become 



3V2 



-a (K) a - b b + 2c ® c) 



De = -p (a(X)a-b(g)b). 
v6 



(24) 



12^/3 



{ (a • r2) a ® a + (b • r2) b b + (c • r2) c (g) c} . 
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FIG. 1: Configuration of the spacecraft constellation in detector's rest frame. 

Note that the expression for detector tensor is higher order in /, compared to the other detector tensors. 

IV. SPHERICAL HARMONIC ANALYSIS OF ANTENNA PATTERN FUNCTION 

The correlation analysis described in Sec. ^1 reveals that the signal modulation induced by detector motion can 
be used to extract the information of the anisotropy of the gravitational-wave background. One important remark 
is that the map-making capability crucially depends on the antenna pattern and/or the detector response function 
in detector's rest frame. We then wish to clarify the relationship between the antenna pattern functions and the 
directional sensitivity to the gravitational-wave backgrounds. To investigate this issue, the spherical harmonic analysis 
of the antenna pattern function is employed and the general rules for multipole coefficients are derived based on the 
geometric properties of the antenna pattern. 

A. Angular power of antenna pattern function 

Similar to the expression ©, antenna pattern function T defined in detector's rest frame is written as 

Tij^L ^) = e^f""-^-'--''^ f)Ffi^^ /) Qif) 

A= + .x 

F/(f2,/)=D,(0,/):e^(0). (25) 
The multipole coeSicient aim for is 

aUf) = r dd r deb sin eY;„,{e, ^) T{f, e, <t>) (26) 
Jo Jo 

with / being defined in (|16|) . We are primarily concerned with how the directional sensitivity depends on the choice of 
the interferometric variables. For this purpose, the optimal filter Q{f) appearing in the antenna pattern function © 
is ignored hereafter. Using the fact that the relations J-*{f, fl) = J-{—f, CI) = J-{f, — O) always hold^, one obtains 

a.,™(/,<) - {-iral_^{-f,t) = {-ir+'al_m{f,t). (27) 
where we used r£_,„(^^,0) = (-l)™r/„(^, </)). 



^ The last equality comes directly from the properties of the respective detector tensors, and thus it holds only among the same types of 
TDI variables, e.g., Ts^Sj ^'^'d TmiM-,- 
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Here and in what follows, we consider the detector configuration in a specific coordinate system shown in Fig. ^to 
calculate the multipolc coefficients. Unless otherwise stated, the unit vectors a, b and c are specifically chosen as: 

V3 1 , V3 1 

a=-— x+-y, b = -y, c = — x + -y. (28) 

While the explicit form of the multipolc coefficients aim depends on the coordinate system (|28|l . a convenient quantity 
invariant under a Euler rotation of the coordinate system can be exploited: 

2£TT ^ (29) 

which characterizes the contribution of ^-th moment to the antenna pattern function. Thus, under the rigid adiabatic 
approximation, the angular power of the antenna pattern in the ecliptic frame is equivalent to that in detector's rest 
frame: 

Ei«^™(/)i' = Ei«^™(/'^)i'- (30) 

rn m 

We use this coordinate invariant quantity to quantify the directional sensitivity of the antenna pattern. 



B. Low- frequency limit 



Consider first the simplest case, / — f / f* ^ 1. In this case, only the £ = 0, 2 and 4 moments for antenna pattern 
function become non-vanishing. This is mostly general except for the fully symmetrized signals such as T-variable. 

In the low- frequency limit, the detector response functions derived in previous section generally becomes of the 
form (see Eqs.ESlEl): 



D 



Ta(a( 



rb(b. 



Tc(c. 



(31) 



except for the T-variable. While the factors T^^b.c may be written as functions of frequency, they do not depend 
on the directional angle fl. Thus, the detector tensor D loses the directional dependence. This means that the 
directional dependence of the antenna pattern function !F arises only through the polarization tensor, e*(n). Since 
the polarization tensor is described by the quadrature of direction vectors u and v (Eq.^), the antenna pattern can 
be generally written as the forth order polynomials of (cos ^, sin 0) and {cos (f>, sine/)). For example, the low- frequency 
limit of the self-correlated signal for Michelson and Sagnac interferometries extracted from the spacecraft 1 is 



/ ^ Ml Ml 



9 



/«1 



-(3 + cos 40) cos^ 9 + {l + cos"* 9) cos^ 0sin^ 



(32) 



Applying the spherical harmonic expansion (|26|l . non- vanishing components of the multipolc coefficients become 



Ml All 

^00 



3V^ 



Ml Ml 



Ml Ml 

40 



The coordinate-free quantity <7i{f) is thus evaluated as 

MiMi 



Ml All 



6v^ 
35 ' 



70 ' 



All Ml 



All Ml 

^4, ±4 



35 



(33) 



for self-correlated Michelson signals. The angular power of self-correlated Sagnac signals are related to that of the 
Michelson signals by cr^^^^ = (4/9)/^ cr^^^*^\ As for the cross-correlated signal extracted from 1 and 2, the antenna 
pattern function is explicitly written as 



/ 3~ Mil 



- ^-^SiSa 



/«1 



-3 
512 



35 + 28 cos 26* cos 46* - 8 (cos 4(/i sin ^) sin'' 9 



f + 0{f). 



(34) 



We then obtain the non-vanishing multipolc coefRcients: 



^00 



10 ' 



'■20 



TT 



-3-1/-, a 



•^40 



140' 



^4, ±4 



lT«\/3 Fk' 
4 V 70' 
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Correspondingly, the invariant quantity ai becomes 



M1M2 _ Sy'tF mjM2 _ _miM2 _ 1 / 477r 



(Tn = , (To ^ ' — , a. — -1/ 35 

° 10 ^ 35 ^ 140 V 3 

for Michelson signals and the same relation cr^^^ ~ (4/9)/^ cr"^"^ holds for cross-correlated Sagnac signals. 

The above examples show that the multipole coefficients ai.„i higher than £ = 4 vanish at the leading order in /. 
Further, the lower multipole moments ^ — \ and 3 also vanish because the antenna pattern ^(Jl) is even function of 
{0,(j)). This is irrespective of the choice of the coordinate system. Indeed, the same properties hold for the optimal 
combinations A, E and T, since these are constructed from the linear combination of Sagnac variables. 

In Appendix ^ the angular power of the optimal combinations are calculated analytically up to the forth order 
in /. It is shown that the lowest order calculation for self-correlated signal cr^* exactly coincides with that for cr™, 
which is related to the self-correlated Michelson signal as af^ — af^ — (2/3)/^cr"^"\ On the other hand, the lowest 
order contribution to the self-correlated signal for T-variable becomes vanishing due to the symmetric combination 
of Sagnac variables. That is, the higher-order contribution of 0(/^) terms becomes dominant in the angular power 
of antenna pattern. While the resultant non-vanishing components for aj"^ are £ = 0, 4 and 6, it turns out that the 
dominant noise contribution for T-variable appears as 0{f^) (see Sec IV C|l . Therefore, the self-correlated signal for 
T-variable is dominated by the instrumental noise in the low-frequency limit and the gravitational-wave signal could 
not be resolved [Ullllla. 

C. Parity invariance in antenna pattern 

Apart from the low-frequency limit as simple limiting approximation, the analytical calculation of ai becomes 
intractable and the perturbative expansion for / generally breaks down. In contrast to the ground-based detectors, 
the difficulty in the space interferometers arises from the transfer function T that appears in equation (|14|l . which 
explicitly exhibits both the frequency and the angular dependences. Thus, to evaluate the directional sensitivity, 
the numerical treatment is required for spherical harmonic analysis. Nevertheless, some important properties in the 
multipole coefficients of an antenna pattern can be still drawn analytically, from symmetric properties of an antenna 
pattern, which is closely-linked with the parity invariance of a detector response. 

Let us introduce two operators, Q : — > tt — and ¥ : (p ^ (p + tt. A composite operator PQ represents the 
parity transformation. The transformation properties of the spherical harmonics Y^ra{0,4>) for the operators are 
fYim = (— l)™y«m and QYim = (— l)^'''™y«m, and so the parity transformation property is PQYfm = {—lYYem- If 
the response function is parity invariant, Eq. (|25|l can be written as 

d9 sin f d(f)Y;^[T+{-lY PQT]. (36) 
Jo 

Then we see that aem becomes vanishing for all odd multipoles {£ — odd) if a response function is parity invariant. 
A similar argument holds for the respective operators P and Q. If a detector response function is invariant under 
the operation P, then aim vanishes for m =odd moment. For response function invariant under the operation Q, the 
multipole coefficients aim vanish when £ + 111 = odd. Here, we summarize these results^: 

aim — for PJF==J^, m = odd, (37a) 

flfcj = for Q = J^, £ + m = odd, (37b) 
airn = for Q¥T = T, £ = odd. (37c) 

Note that while the symmetric property of the antenna pattern itself is a coordinate-free notion, the results presented 
in equations H37a() and (j37b|l depend on a choice of the coordinate system, since the mode m can be mixed by the 
Euler rotation. On the other hand, the property l|37cf) that only depends on £ preserves under the Euler rotation. 



^ There is another interesting property of the multipole coefficient. Eq. I27i tells us that the multipole coefficients of £ = even (odd) 
modes are even (odd) functions of /: 

ai,n{-f) = {-^fai^U)- 

We will see this property explicitly through the low-frequency approximation in the following sections. 
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Keeping this remark in mind, based on the specific configuration and the coordinate system shown in Fig^ and 
Eq. H28|l . several useful formulae related to the parity transformation are derived in Appendix Using these results, 
one finds that the antenna pattern functions for the self-correlated signals of Michelson, Sagnac and the optimal TDI 
variables are invariant under the following three transformations: 

QP^//-^//, ^Tn = Tn. ^Tn^Tu- (J = M„ S„ A, E, T) (38) 

Thus, the multipole moments of antenna pattern functions for self-correlated signals follow the rule (|37|) . This is 
generally true in detector's rest frame under both the static and the equal-arm length limit. The antenna pattern 
function for the cross-correlated Sagnac signals obeys 

^VTij^\Ti.t\\ VTij^\Tij\\ ^Tij^Tij (/,J = S,). (39) 

Several remarks concerning the optimal combinations are in order at this point. First recall that the antenna pattern 
functions constructed from the signals A, E and T can be represented by a sum of the self-correlated and the cross- 
correlated Sagnac signals [see Eas. l|A2| ) (|A3|) ]. Thus, owing to the fact l|39() and the property [jFg.s ]* = s^j it can be 
shown that the property H38I) holds for the self-correlated signals J-aa, -^ee and J^tt, while the cross-correlated signals 
•?^AE, -^^AT and JFet only possess the symmetry, = T . Hence, the cross-correlated data may generally contain the 
I = odd moments. Note, however, that in the low-frequency limit, the appreciable multipoles are the ^ = 0, 2 and 4 
modes. This readily implies that the contribution of the i = odd modes becomes significant at / = ///* ~ 1, which 
will be explicitly shown in next Section [see angular power in Fig. |31and effective sensitivity curves in Fig.0]. 



D. Geometric relation between optimal combinations of TDIs 



In a specific case with signals constructed from the optimal combinations (A, E,T), a further important property 
is obtained combining with the geometric relationship among the three spacecrafts. 

Let us start with the property of the Wigner D matrices ()l()|l . For specific choice of the angles i9 = and = tt, 
the Wigner D matrix becomes (i^„(0) = bnm and d^^(7r) = (— l)^+™(5„,_m, respectively js^- Thus a rotation i/; with 
"i? = <y9 = transforms the coefficients ag^a to as 

= e-™'^af„,. (,? = = 0) (40) 

On the other hand, a rotation -i/; followed by the rotation -d = if = -k transforms the coefficients as 

4™ = (-ire™'^Km]*. (z?-^-^) (41) 

From the spacecraft constellation shown in Fig. ^ the antenna pattern functions of the self-correlated Sagnac 
signals T^^^^ and are related to T^^^^ by the Euler rotation angles i/' = 27r/3 and ?/' = 47r/3 — ip = 0), 

respectively. Similarly, all multipole coefficients of cross-correlated Sagnac signals are related to a^,^j^ as indicated by 
(|40|l and H41|l . These relations are summarized as follows: 

"£m ~ ^ "£m ' "'Ira ~ ^ "'£771 ' 

„SiS3 / 1 \mr^SiS2i=^ ^^2^1 / 1 \m^ — 4im(5r^SiS2]* ^^3^2 / 1 Am „ — 22m(5 r^Si S2] ^ 

where 5 = 27r/3, and we have used l|27|l . This means that the antenna patterns for all the possible combinations of 

^L'' and ag^l^ 



(A, i?, T) can be represented by a sum of the primary multipole moments of the Sagnac signals, a^^^ and a^^-^- In 



this sense, the optimal combinations of TDIs are not strictly independent. 

Based on the relations 142|) . with a help of the expressions (|A2p . the multipole moments for self-correlated signals 
J^AA and J^EE can be rewritten with 

= e-"^ a^^::^ - a^^^ + (-1)" [a,^^=]*} , 

af- = D(^^ e-^"^' - Di?,^ {e-'^"' + (-1)" [a^'l*} - (43) 

where the coefficients Cm and Dm become 

CW=cos(m<5), = i = i{2-Fcos(m<5)}, Z?L'^ = i {4cos(mJ) - 1} . 
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TABLE I: Symmetric properties of antenna pattern function 



Combination of variables 



Condition 



Properties of (/) 



AU 
All 

(A,A), (E,E) 

(A,E) 
(A,T), (E,T) 



low-frequency limit (/ <^ 1) 
self-correlation 



for £ / 0, 2, 4 
for ^£ =odd 

= 

for ^ = 0, ft 
for £ = 



The details of the proof are presented in Appendix lO 



Now recall from the properties l|37(l and (|H8(I that the non-vanishing components of the multipole moments (|43|l are 
the i = even and m = even modes. Then the comparison between the coefficients Cm and Dm leads to the relation 
«£m = <m for TO = 0, ±6, ±12, ±18, • • • and a^* = -a™ for m = ±2, ±4, ±8, • • • . Thus, one finds 



ar=ar. (44) 

The similar identity also holds for the cross-correlated signals Tat and JF^t- Applying the relation (|42|l to the 
expressions HA3|I . one obtains 



2jv/2 _,„,,, ^ fmS 
3 

= tan j aZ^. (45) 

Thus, both of the multipole moments af^^ and aj^ become vanishing when to = 0, ±3, ±6, • • • . Further, for all the 
non-vanishing components, the absolute value of the pre-factor becomes unity. This immediately yields the relation: 

af = oY- (46) 

Note that while the relations H44I) and 1)46(1 are derived in a specific choice of the coordinate system (|28|l . the final 
results do not depend on the coordinates. 

Finally, we note a quite remarkable fact for the cross-correlated signals, T t^^^ Tj^^- and T^^^. It is shown in Appendix 
121 that the multipole moments £ = Q and 1 for the antenna pattern are exactly zero, while the monopole mode 
[l = 0) vanishes for the cross-correlated signals T^t and Tet, over the whole frequency range: 

o'o(/) = 0, o-i(/) = for AE-correlation, 

CTo(/) = for AT-,ET-correlation. (47) 

Here, the important properties of the antenna pattern functions derived from the parity invariance and geometric 
argument are summarized in Table 



V. ANGULAR POWER AND DIRECTIONAL SENSITIVITY OF SPACE INTERFEROMETER 



While several important properties for directional sensitivity of the space interferometer were found in previous 
section, it remains still unclear how the multipole moments of the antenna pattern functions quantitatively depend 
on the frequency and the combinations of data streams. In this section, based on the previous remarks, the spherical 
harmonic analysis of the antenna pattern function is carried out analytically and numerically in specific choices of the 
data combinations. For a relevant range of the frequencies beyond the low-frequency approximation, the directional 
sensitivity to the antenna pattern is estimated in the LISA case, taking into account the instrumental noises. 



A. Toy model example 

As noted in Sec lIV Cl the frequency and angular dependences of the transfer function T in equation H14() make 
it difficult to treat the spherical harmonic analysis of the antenna pattern. If we set T = 1, however, the spherical 
harmonic expansion of the antenna pattern can be analytically evaluated, the results of which are compared with the 
realistic cases without invoking the assumption T = 1. 
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For computational purpose in this subsection, we set the directional unit vectors for three spacecrafts as: 

a = z, b=— y--z, c = — 2""^ ^ 2^' ^ ' 

and consider the antenna pattern for Michelson signal. Table suggests that a number of non- vanishing multipole 
moments is severely restricted in the case of the self-correlated signals, since the assumption T = 1 roughly corresponds 
to the low- frequency limit. An interesting case is therefore to take the cross-correlation between the signals extracted 
from the vertices 1 and 2. In this case, the response function at the rest frame becomes 

A=+,x [ = i[b®b-a(»a] : e^(ri) 

With the specific choice of the coordinate system (|48f) . the explicit expression for the antenna pattern becomes 

^12{IM) = -—\icos^<l) + 2sin^ecos^<i) + ism^e 
64 

+ cos^ e sir? (/) {-2 + 5 cos 26* + 2 cos 20 + 3 cos^ 9 sin^ 0}] e"'^ ™" (49) 
Note that the function H49|l possesses the following symmetry: 

QP.F12 = [Tl2\* , P^12 = ^12, Q^12 = [^12]* . (50) 

This indicates that the antenna pattern is sensitive to both the even and the odd modes, while the multipole moments 
aim with TO — odd become vanishing. Since the relation H27I) always holds, it is sufficient to treat the case for m > 
only. Substituting the explicit expression into the definition (f^ . the integral over </> is first performed. Writing 
cos by u, we have 



where the function gm{u) can be expressed as polynomial series as 

g„(^)=^c(™'^)u2A^(l-u2)H/2. (52) 
Af=0 

The coefficients c^™'^^ are the numerical constants, which are summarized in Appendix IdI Note that the function 
are non-vanishing only for to = 0, ±2, ±4, indicating that the non-vanishing components of aim are obtained only 
when TO = 0, ±2, ±4. From l|51(l and (|52|l . the remaining integrals become of the form: 



This integral is analytically performed according to Ref. |23|- Using the formula for Legendre polynomials, P™{u) 
(— 1)™(1 — v?)"^/'^{d/du)"^Pii{u) for to > 0, repeating the integration by parts yields 



rN 



rl™ r 



f_i du 

2(m+JV) 



E 3^'"^^/) / duPi{u)u^e-^f\ (54) 

The quantities gs™'^\f) are the polynomial function of / up to the forth order and are listed in Appendix IdI The 
integral in the last line is expressed in terms of the spherical Bessel function jg: 

duPi{u)u'e~'f'^ = i'-^ I du Pi{u)e-'f'^ 
-1 df'^ J -I 

= 2(-l)^ Z+^-C- nU). (55) 
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Thus, substituting the results (|52() . H54() and H55|l into H51() . one finally obtains the analytic expression for multipole 
coefficients: 

/————— — 2 2(rn+N) 

While the above expression is the outcome based on the coordinate (|48|l . the invariant quantity (Ji{f) can be evaluated 
from (|56|l . which is depicted in Figure |2] as a function of £ and /. Also using (|56|l . the non- vanishing components of 
cr^ in the low-frequency limit are explicitly calculated as 

^ _ /2 ^ \^ ? _ 17%/2^ .3 ^ 3V¥ (V30-6)V¥ ~2 

10 28 ' 14 2520 ' 35 960 ^ ' 

1 /19^ ^ 43 t?, 1 /47^ 19 /T9^ y2 

^^~42VlO'^ 1848V 190''' ~ 14oV 3 264V 141 ^ 

up to the third order in /. The leading order terms in ui rigorously match the results in equation H35|l . 

Figure 12 shows that the higher multipole moments appear as increasing the frequency, and an oscillatory behavior 
is found in the frequency domain / ^ 1, which are also indicated from the low- frequency expansion in equation |S7J. 
From the analytic expression (|56|l . we readily see that the quantity higher than ^ > 4 scale as 0{f^~'^) in the 
low-frequency limit and asymptotically behaves as ai oc /^^ in the high-frequency limit. The resultant angular power 
depicted in Figure |2] implies that the resolution of anisotropy in the stochastic background of gravitational waves can 
reach ^ < 10 ~ 15 for a relevant frequency range < / < 10. This result is comparable to the angular resolution of 
gravitational- wave background measured from the ground-based detectors [23 . |28| , since the assumption neglecting 
the frequency and the directional dependences of transfer function T can be validated for the response function of 
Fabry-Perot interferometer. 
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FIG. 2: Angular power o"^(/) of the antenna pattern function for the toy model as a function of £ and /. 



B. Influence of transfer function 



We now calculate the angular power of the antenna pattern fully taking into account the frequency and the angular 
dependences of the transfer function T. The results are then compared with the toy model example. For this purpose, 
the spherical harmonic expansion of the antenna pattern function is numerically carried out using the SPHEREPACK 
3.1 package j3"§ |. 



14 



FigureOlshows the angular power of the antenna pattern function for the self-correlated Michelson signals J-^^m^ (Isft) 
and the cross-correlated Michelson signals ^miM2 (right). Relaxation of the assumption T — const., i.e., low-frequency 
approximation, leads to the non-vanishing components for even modes with £ > 6. However, the resultant higher 
multipole moments turn out to be highly suppressed. While the sensitivity to the higher multipole moments is 
slightly improved in the case of the cross-correlated Michelson signals, comparing it with Fig. |2 reveals that the 
frequency dependence of the transfer function T significantly reduces the angular power in both the lower and the 
higher multipole moments. The numerical evaluation of spherical harmonic expansion implies that the non-vanishing 
components of the angular power asymptotically decrease as ai oc /^^ in the high frequency region, even faster than 
that of the toy model example. 

The behaviors of the angular power are qualitatively similar to the case adopting the Sagnac variables that cancel 
the laser frequency noise (Fig0}. Apart from the low-frequency limit, where the antenna pattern function for Sagnac 
signals behaves as ^s^s^ ~ j [see Eqs. and ^^], the angular power is highly suppressed at the frequency / > 1 
even in the relatively lower multipole moments ^ 6. Thus, the directional sensitivity of the space interferometer to a 
stochastic background of gravitational waves is severely limited by the frequency dependence of the transfer function. 
This fact is irrelevant to the choice of the interferometric variables. 




FIG. 3: Angular power creif) of the antenna pattern function for the self-correlated Michelson signals J-MiMi{left) ^-nd the 
cross- correlated Michelson signals, ^miM2 (right). 



C. Directional sensitivity for optimal combinations of TDI variables 

To elucidate a more quantitative aspect of the directional sensitivity to the gravitational-wave background, it will 
need to take into account effects of detector noises. To investigate this, rather than using the Michelson and Sagnac 
signals, a set of optimal TDIs (A, E, T) free from the noise correlations should be applied to the correlation analysis 
of the gravitational-wave signals. 

Figure plots the quantity ai for various combinations of the optimal TDIs. Note that the angular power of 
the antenna pattern function JFee (-^et) coincides with that obtained from J^aa (-^at), although the sky patterns 
themselves differ from each other. For the self-correlated signals, the amplitude of JFaa is quite similar to that of 
the self-correlated Sagnac signals ^sjSi : while the low- frequency part of the angular power for JF^t is highly suppressed, 
which can be deduced from the low-frequency approximation presented in Appendix^ As for the cross-correlation 
signals, the monopole and the dipole moments for the antenna pattern function J^^b are exactly canceled and the 
monopole moment for further vanishes (Tableland Appendix [2)). Apart from these facts, the magnitude at 
frequency / > 1 shows a rich structure with many peaks, indicating that the directional sensitivity could be improved 



15 





FIG. 4: Angular power o-e{f) of the antenna pattern function for the self-correlated Sagnac signals ^s^s^ {left) and the 
cross-correlated Sagnac signals Ts^s^i^^ght) ■ 

at / > 1. As shown in Fig. the angular power of the cross-correlated signals is one order of magnitude smaller 
than that of the self-correlated signals, however, this does not directly imply that the self-correlated signals are more 
sensitive to an anisotropy of a gravitational-wave background. 



Based on these results, let us now quantify the directional sensitivity of the antenna pattern function. Assuming 
that the laser frequency noise can be either canceled or sufficiently reduced by the TDI technique with the use of the 
recently proposed laser self-locking scheme 39, 40], the dominant noise contributions to detector's output would be 
the acceleration noise and the shot noise. According to Ref. 33], the noise spectral densities for optimal TDIs are 
calculated as (see also 



^r(/) = ^r(/)=sin'(/72) {S (2 + C0S/) 5shot(/) + 16 (3 + 2 cos / + cos 2/) 5accci(/)}, 
+ 2cos/) {5shot(/) +4sin2(/72) 5 

acccl (/)}■ 



sru) 



1 



(58) 



Note that the cross-correlated noise spectra are exactly canceled. Here we specifically adopt the noise functions for 
the LISA detector: S'shot(/) = 4.84-42Hz-i and S^ccdlf) = 2.31-*°(mHz//)4Hz-i|33|. We then define the effective 
sensitivity for the multipole moment £, h^J'^{f), which characterizes the rms amplitude of the noise-to- angular power 
ratio: 



ISM) 



(59) 



for self-correlation signals. Setting i = 0, the above definition recovers the usual meaning of sensitivity curve. Thus, 
the quantity h^^^(f) may be regarded as the effective power of £-th moment relative to the monopole moment as a 
reference sensitivity. For the cross-correlated signals, on the other hand, the absence of noise correlation implies that 
the signal-to-noise ratio can be improved by optimally filtering the cross-correlated signals. The resultant form of 
the signal-to- noise ratio shows the explicit dependence on the observation time T ^3]- According to Ref. j2^], the 
effective sensitivity for cross-correlated signals may be written as: 



(ta/) 



1/4 



SnAif) Sr^M) 



1/4 



(60) 



where A/ denotes the frequency resolution for actual output data. 
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FIG. 5: Angular power o-e(f) of the antenna pattern functions for the optimal TDI variables. The top panels show the 
magnitude ae in the self-correlated cases TAA{left) and J^tt (right), while the bottom panels represent the results obtained from 
the cross-correlated signals J^AE{left) and J-AT{right). Note that the angular powers of antenna pattern JTee and J-^t coincide 
with those of J^aa and J-at- 

Figure shows the effective sensitivity curves for the self-correlated and cross-correlated optimal TDIs as functions 
of / = ///*• In plotting these curves, we used the characteristic frequency /* ~ lOmHz for the LISA detector. The 
different lines in each panel indicate the effective strain sensitivity for each multipole moment. Clearly, the directional 
sensitivity to a gravitational-wave background is not so good in the case of the self-correlated signals. As anticipated 
from Figure O and the noise spectra H58|) . the effective sensitivity in the low- frequency limit scales as h^fg cx for 

£ = 0, 2 and 4 of AA-correlation and h^^g oc /^'^ for £ = 0, A and 6 of TT-correlation. At the frequency around 
the characteristic frequency the directional sensitivities may reach at a maximal level and the higher multipole 
moment can be observed in both AA- and TT-correlations, however, the detectable multipole moments are still limited 
to £ =even mode with £ < 6 for the sensitivity h'fg ~ 10^^* Hz^^/^. 

The situation might be improved if we consider the cross-correlation signals. In bottom panels of Figure |S1 the 
observation time of T = 1 year and the frequency resolution with interval A/ = //lO are assumed. In this case, the 
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sensitivity reaches h^g ~ 5 x 10~'^^llz~^/'^ in both AE- and AT-correlations (and also the ET-correlation) , and the 
detectable multipole moments become, say, ^ < 10 or even higher multipole moments in both £ =odd and even modes. 
At the frequency ///* ~ 3, the effective sensitivity for the higher multipole moments becomes comparable to that for 
the lower multipole and shows a complicated oscillatory behavior. Although the antenna pattern for cross-correlation 
signals is completely insensitive to the £ — mode, improvement of the sensitivity is noticeable, which might be useful 
to distinguish between the gravitational-wave backgrounds from Galactic origin and those from extragalactic origin. 




FIG. 6: Effective sensitivity curves defined in H59|l and H6()^ for the self-correlated and cross-correlated optimal TDI variables: 
AA-correlation(top-/e/f); TT-correlation(top-ri(;/if); AE-correlation(6oitom-Ze/f); AT-correlation(6ottom-rjg/i£). In cases of the 
cross-correlation signals depicted in bottom panels, the observation time of T = 1 year and the frequency resolution A/ = //lO 
are assumed. The characteristic frequency is /« ~ lOmHz. 



VI. CONCLUSION & DISCUSSIONS 



In this paper, we discussed the directional sensitivity to the anisotropy of gravitational-wave background observed 
via space-based gravitational-wave detector. While the detection of anisotropic gravitational-wave background could 
be achieved utilizing the modulated signals of cross-correlated data induced by the detector motion, the directional 
sensitivity and the angular resolution crucially depend on the antenna pattern function and/or the detector response 
in detector's rest frame. In contrast to the groundbased detector, the space interferometer with long baselines gives 
a rather complicated response to the gravitational-wave signals. 

We have performed the spherical harmonic analysis of antenna pattern function for space interferometer and studied 
the general features of antenna pattern sensitivity beyond the low-frequency approximation. We have shown that 
the sensitivity to the multipole moments of an anisotropic gravitational-wave background is generally restricted by 
the geometry of the detector configuration and symmetries of the data combinations (see Table Sec. IIVCI and 
Sec. II V l51l . The numerical analysis of the antenna pattern functions reveals that the angular power of the detector 
response increases with frequency and shows the complicated structures. To characterize the directional sensitivity, we 
introduced the effective sensitivity /i^.^ (/) for each multipole moment and evaluated it in the case of the LISA detector 
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specifically. Using the cross-correlated data of optimal TDIs, i.e., AE-, AT- and ET-correlations, we found that the 
detectable multipole moments with effective sensitivity /i^^ 10~^°Hz~-'^/^ may reach £ = 8 ~ 10 or even higher 
multipoles at / ^ /* = 10 mHz, which would be useful to discriminate between the gravitational-wave backgrounds 
of Galactic origin and those of the extragalactic and/or the cosmological origins, recently discussed by several authors 
(e.g., Illlllfeo]). 

Although the improvement of the directional sensitivity beyond the low-frequency approximation is remarkable, 
the sensitivity of the space interferometer is still worse than the one achieved by the cosmic microwave background 
experiments, like the COBE (cosmic background explorer) and WMAP (Wilkinson microwave anisotropy probe). The 
one main reason is that the wavelength of the gravitational waves to which the space interferometer is sensitive is 
comparable to or longer than the arm length of the detector. Because of this, the response to the gravitational-wave 
background becomes simpler and most of the directional information is lost, as seen in Sec UVBl The directional 
sensitivity can be improved as increasing the frequency, however, the sensitivity beyond the characteristic frequency 
/* is limited by the instrumental noises. Another important aspect is that the phases of the gravitational-wave 
backgrounds are, in nature, random. Thus, the information of phase modulation induced by the detector motion 
cannot be used. This is marked contrast to the signals emitted from the point source, in which the angular resolution 
can reach at a level of a square degree or even better than that [3, S IHi HE uM ■ 

Further notice the important issues concerning the map-making capability of the gravitational-wave backgrounds. 
As discussed in Sec^ provided the time series data, the task is to solve the linear system (lllll under a prior knowledge 
of the antenna pattern functions for the space interferometer. The crucial remark is that the antenna pattern functions 
for the cross-correlation signals taken from the optimal TDIs (A, E, T) are not independent (see Sec lIVD]l . This fact 
implies that Eq. Hll|) constructed from the three cross-correlation data (AE, AT, ET) is generally degenerate. Thus, 
the deconvolution problem given in 1)11(1 would not be solved rigorously. Rather, one must seek a best-fit solution 
oi pf^{f) from Hll|) under assuming a specific functional form of the luminosity distribution p^„(/). The analysis 
concerning this issue is now in progress and will be presented elsewhere. 
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APPENDIX A: SPHERICAL HARMONIC EXPANSION IN THE LOW-FREQUENCY 

APPROXIMATION 



In this appendix, employing the perturbative approach based on the low- frequency approximation / = ///* ^ 1, 
the spherical harmonic expansions for several antenna patterns are presented, which partially verify the properties 
summarized in Tableland Eq. ()37fl . 



1. Sagnac interferometers 

The angular powers of the self-correlated Sagnac signal J-s-s^ and the cross-correlated Sagnac signal Ts-sj in the 
low- frequency approximation are summarized as follows: 
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2. Optimal combinations of time-delay interferometry 

From Eq. H22() the antenna pattern functions of the self-correlated optimal TDIs can be written down in terms of 
the antenna pattern for the Sagnac signals: 

^ AA — 2 {"^ ^33 + -^ll ^ (13)} 1 

•^EE — 77 {-^ 11 + 4^ 22 + 33 + ^ (1,3) (12) 2^ (2,3)} J 
1 

•^TT = 2 {-^ 11 + -^22 + ^33 + (1,3) + ^ (1,2) + (2,3)} ■ (-^2) 

where the round brackets are an abbreviation for the symmetrization, ^(1,3) = + ^31 for instance, and J-^^ stands 
for ^s^sj- The antenna pattern functions of the cross-correlated signals, J-ae, ^at and J-'eti are 

^ AE 2\/3 ^ 11 ^" *^3i ^ 13 2^ 12 2^ 32} 5 

^ AT ~ ^ {*^ 33 11 ^ 31 13 ^" *^ 32 *^ 12} 7 

V6 

^ET - ^ {^11 - 2^22 + ^33 + ^12 - 2^21 + ^32 " 2^23 + ^(1,3)} • (A3) 

Under the configuration in a specific coordinate H28() , the low- frequency approximation of the multipole coefficients 
of the self-correlated optimal variables are given as follows: 



^ee: aoo-^/, a2o--^-/, «« - ^/ , «44 - (^g ^ ^ + ^ ^ ^ j / , 

Tt.: «oo-^/, «40--y^/, "^^-^i^TH^- ^^^^ 

TO < modes are given by the relation (|27|l . The rule H37|) strictly restricts the appearance of multipole moments, 
and of course the above multipole moments follow the rule (|37() . For the cross-correlation of two data streams, one 
would expect that i = odd modes appear even in the low-frequency limit. However it is not the case. A non- vanishing 
multipole moment is given to order 0{f^) by 

•^AE : 044 = ( -i= - -4= ) V^P- (A5) 



,2V210 6V70 

The i = odd modes appear in the next order 0{f^) in some multipole moments that satisfy £ + m — even: 



:Fae: a33 = 3^^|/^ «53 = ^J^/^ (A6) 



The angular powers up to 0{f ) are summarized as 



20F^^2 211^^4 _4v^-2 16\^?4 _2v^-2 4511^,4 _ /1829 .4 



^aa,^ee: <to- ^ / 1260-^' 35 315-^' '^""105'^ 498960'''' '^'^ ~ 72072 V 15^ 



: o-o = —r, (74 = — — ere = - — ^^f* (A7) 

" 504-' ' 4752-' ' ^ 144144 ^ ' 



and 



T n n ^P 

■Tae ■ Co = 0, (71= U, (72 = p, (73 = =, 

126V3 9\/30 

"^"^^"WV7^' ^^~WlTn' ^'"WlTB- ^^^^ 
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APPENDIX B: PARITY TRANSFORMATION 

Here, we summarize some formulae related to the parity transformation, which are used in Sec lIV CJl In general, 
parity of the polarization tensor e+^^ depends on the choice of the coordinate basis. Our choice of the basis vectors 
are those defined in Q and I0J just simply replacing the variables dE,4'E with d,(f> in detector's rest frame. In the 
following, the vector d stands for the unit vectors a, b, c. We then obtain 

P(d • rj) = -1 (d • O), 

P(d(g)d) : = +1 (d(8)d) : e-^ {A^+,x) (Bl) 



for the operator P and 

Q (d (g) d) : e+ 
Q(d(g)d) : e"" 

for the operator Q. As for the composite operation QP, 

QP(d-rj) 

QP(d(g)d) : e+ 
QP(d® d) : e"" 



= +1 (d-fi), 

= +1 (d(g)d) : e+, 

= -1 (d® d) : e"" (B2) 



= -1 (d-n), 

= +1 (d(Xi d) : e+, 

= -1 (d® d) : e^. (B3) 



which is identical to the parity transformation, one has 



APPENDIX C: ON CANCELLATION OF MONOPOLE AND DIPOLE MOMENTS IN ANTENNA 
PATTERN FUNCTION FOR CROSS-CORRELATED OPTIMAL TDIS 

In this appendix, we will prove that the antenna pattern function for the cross-correlated optimal TDIs has the 
following symmetric properties: 

o'o(/) = 0, cTiif) = for AE-correlation, 

cro(/) = ' for AT-,ET-correlation, (CI) 

which are intimately related to the geometric properties of both the detector configuration and the response function. 
As we have explained, the antenna pattern functions J-'ae, ^at and J-^t are written in terms of the antenna pattern 
functions for the Sagnac signals [see Eq. (|A3p ]. The expressions readily imply that the multipole moments for the 
antenna pattern functions, af^, af^ and a^^^, are also obtained from the sum of the cross-correlated Sagnac signals, 
^em' [see (^Sj, for example]. 

Let us first consider the monopole moment aoo- Since the monopole moment is obtained through the all-sky 
average of the antenna pattern function, it is, by construction, invariant under both the Euler rotation and the parity 
transformation of the coordinate system. This indicates that the monopole moments for various combinations of the 
Sagnac signals a^^/ are degenerate and there are only two independent variables, that is, 

aoo''' = "m^ = Ooo''': "m' = ^00''' 7^ 0- (C2) 



Substituting this into the spherical harmonic expansion of the antenna pattern functions (|A3|) , we immediately see that 
the monopole component of cross-correlated optimal signals exactly vanishes, i.e., ago — Oqq — ago = 0- Accordingly, 
the monopole moments of angular power, ctq"^, ctq^ and CTo^, become vanishing. 

Next focus on the dipole moment of AE correlation. In a specific choice of the coordinate system (|28ll . all the 
components in the dipole moment vanish due to Eq. (|37|l for the self-correlated Sagnac signals, i.e., a^'^^ = 0. Also, 

the dipole moment with m = becomes zero for the cross-correlated signals, i.e., a^o^^ = 0. Further, the relation H27|) 
implies a^lj^ = [a*f]*. Collecting these facts, the dipole moment of angular power af^ can be written as 

I |a-| ^ ^ \alf^ - ali^^ + 2alf^ - 2a^f ^ | . (C3) 



It is thus sufficient to consider the dipole moment with m = +1 for the cross-correlated Sagnac signals. 
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From Eq. (|42|l , the angular power erf can be solely determined by the quantity afl^^ . If we write al\^^ = r{f) e^^^^^ , 
Eq. (|n3|) becomes 



a^^if) = ^ |2 + e-^^ + (l + 2e-2^) ^-,:2.|^^^^^^_ (C4) 

To determine the phase factor 9{f) or amplitude we recall the fact that the dipole moment of antenna pattern 
function J^aa vanishes: 

AA ^ r SiSi , SoS-, SiSa SgSi-i 1 / S1S3 , S-,Si\ 

= 0. 

from (|22|) . Using the relations (|42|l . the above equation becomes 

AA _ / -129 _ -i25\ _Q 

which finally yields r{f) = or 9{f) ~ 5 + rnr. Thus, substituting this value into the right-hand-side of equation (jC4|l 
immediately leads to the conclusion that the dipole moment of angular power crf^ is exactly canceled. This means 
that all the dipole components for af^ become zero over the whole frequency range. 

APPENDIX D: COEFFICIENTS IN ANTENNA PATTERN FOR TOY MODEL 

In this appendix, we summarize the coefficients c*-™'^^ and the functions gi"^'^\f) defined in H52I) and (|54|l . 
For the coefficients c^™'^\ the non- vanishing components are 









(0,0) _ 1237r 
256 


„(o,i) _ 1357r (o_2) _ 
128 ' 


2197r 
256 ' 








^(2,0) _ IStt 
128' 


c(2,i) _ 397r 
128' 










512' 






■ the functions gi 




the non-vanishing 


components for < ?7i < 4 and < < 2 become 




= 1, 


9^' - 


-4-/^ 


9^'-'' - 2, 


= 144 + 48/2 + /4, 






5f = 


-4*./, 


5f'°U-576z/-32z/3, 






= 1, 


(2,0) 
52 = 


12 + 2/2^ 


#^'=-24-/2^ 


g!^'°'^ = -1440 - 432/2 _ 


9^' 


= -8«/, 


(2,1) 
53 = 


16*/, 


5^"'°^ = 1920i/ + 96i/3^ 




(0,2) 
54 


= 1, 


(2,0) 
54 = 


-/^ 


5f^^=30H-2/2, 


= 1680 + 720/2 _^ 6/^ 


(2,1) 






-1344i/- 96i/3, 






9^' 


= -/^ 


9^ - 


-336/2 _ 4^4^ 






5^"' 


= 32i/3, 










9^ 













(Dl) 



Note that the other components with m > do not contribute to the calculation of the multipole coefficient aim due 
to the coefficients c("'^) (Eq. ffHT^ ]. 
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